Fractional Generalized Langevin Equation Approach to Single-File Diffusion 
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Fractional generalized Langevin equation with external force is used to model single-file diffusion. 
It is found that for external force that varies with power law the solution for such a fractional 
Langevin equation gives the correct short and long time behavior for the mean square displacement 
of single-file diffusion when appropriate choice of parameters associated with fractional generalized 
Langevin equation are used. By considering some special cases of the fractional generalized Langevin 
equation, a new class of closed analytic expressions for the mean square displacement of single-file 
diffusion can be obtained. The effective Fokker-Planck equation associated with single-file diffusion 
is briefly considered. 
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I. INTRODUCTION 

Recent advances in nanofabrication allow the preparation 
of new types of nanotube materials such as carbon nan- 
otubes, petide nanotubcs, inorganic and organic zeolites, 
etc. The study of the molecular transport processes in- 
side these nanotubes or channels have attracted consider- 
able attention [3, S [H [H, [M H3 ■ The one dimension 
transport process of an assembly of non-passing particles 
(without mutual passage) in confined geometries such as 
narrow pores or nanotubes is known as single-file diffu- 
sion (SFD). In other words, particles undergoing SFD 
maintain the order of their arrangement at all times. 

The main feature of SFD can be characterized by its 
short time and long time limits of its mean square dis- 
placement (MSD) denoted by R 2 . For very short ob- 
servation time, the particles diffuse normally and satisfy 
Fick's law, such that the short time limit of MSD of SFD 
is given by 



R 2 



(x(t) 



2D Q t, t ->• 0, 



;i) 



where D is the diffusion coefficient. In other words, 
for very short diffusion time, the motion is just ordinary 
Brownian motion, which is a Markov process. As for the 
long time limit of the MSD for SFD, one has 



R 2 = 2FVi, t -> oo, 



(2) 



where F is the SFD mobility. Recall that diffusion that 
does not satisfy Fick's law is known as anomalous diffu- 
sion with MSD satisfying R 2 oc t a , a ^ 1. When a > 1 
the diffusion is enhanced and it is called superdiffusion. 
In the case with a < 1 the diffusion is subdued and one 
has subdiffusion. Thus, the long time behavior of SFD 
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belongs to subdiffusion, which is a non-Markovian pro- 
cess, indicating the motion is correlated. 

Consider the simple case of a single molecule in a one- 
dimensional pore, such that the molecule can move in ei- 
ther direction with equal probability via activated jumps 
with step length I. If t is the average time between suc- 
cessive jumps, and is the fractional occupancy, then one 
has 



R 2 = l 2 (l 



T 



0. 



(3) 



which corresponds to normal self diffusion. Fedders [li 
has derived the long time limit of the MSD as 



R 2 = I 



(1-0) 



t — > oo, 



(4) 



which shows that at longer times the motion is strongly 
suppressed by collisions with neighboring particles (see 
also [HI). From © and (JU) one obtains the SFD mobility 

as 



F = I 



2 (1- 



1 



(5) 



About a decade ago, Brandani Q introduced heuristi- 
cally (without derivation) the following analytic expres- 
sion for the MSD in SFD: 



R 2 



V 



(1 - 0)t/r 



(6) 



Expression ([5]) is able to give the correct limiting cases 
(|3]) and (01 . Even until today, (|6|) is still regarded as an 
expression that "comprises both cases with satisfactory 
accuracy" [Ref. [13 , page 334]. 

A similar expression for the MSD of SFD was obtained 
by Lin et al [301 ] based on the following ansatz: 



1 

R 2 



1 



1 



2D a t 2Fy/t' 



(7) 



2 



By solving (JT)) for MSD one gets 



1?- 



2D n t 



1 + DoVi/F' 



(8) 



Just like ([5]), one can obtain from ((SJ) the short time 
and long time limits ([I) and ((2]) respectively. To the 
best of our knowledge, so far there still do not exist any 
derivations for the closed analytic expressions ^ or © 
Although the research on SFD can be traced back to 
1960s, but until now there are not many studies that can 
provide a comprehensive description of the process. One 
notable exception is that of reference [l^ , which has de- 
rived the exact SFD propagator valid for all time scales 
based on the reflection principle of Chandrasekhar Q. 
Other attempts to model SFD include the early statisti- 
cal and probabilistic models QJ, H El 11 [H, \M, M, HI 
to the more recent ones based on fractional dynamics 
[3, EL & HH, HH . The main aims of this paper are two 
folds. First we introduce a new type of fractional gen- 
eralized Langcvin equation with external force to model 
SFD. Our second objective is to derive a new closed ex- 
pression for the MSD of SFD. Despite that the fractional 
generalized Langevin equation under consideration does 
not lead to a closed expression for the MSD of SFD, it is 
still possible to show that it gives the correct short and 
long time limits under some specific conditions. It is pos- 
sible to derive from various special cases of the fractional 
generalized Langevin equation a completely new class of 
closed expressions for the MSD of SFD, which can be 
regarded as alternatives to ^ and ©. Our previous 
work [29[ which though contains a detailed discussion 
of fractional generalized Langevin approach to SFD and 
its realization as a step fractional Brownian motion, it 
does not provide the derivation of closed expressions of 
the MSD for SFD. In addition, the present work also pro- 
vides a discussion of the effective Fokker-Planck equation 
for SFD. 



II. FRACTIONAL GENERALIZED LANGEVIN 
EQUATION WITHOUT EXTERNAL FORCE 

In this section we want to study whether it is possible 
to model SFD using the fractional generalized Langevin 
equation without external force. For the purpose for sub- 
sequent discussion, let us first consider the following frac- 
tional generalized Langevin equation with external force: 



D a v(t)+ j( t -T)v(r)dT = f(t)+at), (9) 
Jo 

v(t) = Dx(t), (10) 

where < a < 1. is an internal Gaussian noise with 
mean zero and covariance 



mas))=c(\t- s \) 



(ii) 



fit) is a time-dependent external force given by 

t- K 



/(*) = a 



< K < 1, 



(12) 



and "f(t) is the memory kernel which will be specified later 
on. The fractional derivative used in © is the Caputo 
fractional derivative, which is defined by [H, HH, HU 

D a g(t) = I m ~ a D m g(t), (13a) 

where the fractional integral I a is defined for a > as 

1 



r (a) Jo 



(t - u) a - l g(u)d Ul 



(13b) 



with m — 1 < a < m, and m is positive integer. 

Note that generalized Langevin equation @ with 
fit ) = has been studied by several authors [111 fl2l . 
I29I [3lj . The solution for the position process from © 
and (fTUj) with initial conditions x(0) = x , u(0) = v Q is 
given by 

x(t) = x +vJ 1 - a Git) +aI 1 ~ K Git) + [ Git~u)£iu)du 



(14) 

where Git) is given by the inverse Laplace transform of 

1 



G(s) = — 
The mean of xit) is given by 
x = (xit)) = x + 



■j(s)s' 



(15) 



r(i - a) J 



it-u)- a Giu)du 



r(i - «) Jo 



(t - u)- K Giu)du (16) 



and its variance is 



ixit) - xY 



duGiu) / dvCiu- v)G(v), (17) 



The MSD of xit) in terms of its mean and variance is 



R 2 



(ixit)-x ) 2 ^ = ix-x ) 2 + a 2 . (18) 



Here we shall consider two cases without external force, 
that is with the constant a = 0; the case with external 
force will be studied in next section. The generalized 
Langevin equation (|9|) includes the following cases: 

Case 1. < a < 1, 7 (i) = Xi5(t), Xi > 0, C(t) = cS(i), 
and c > 

© becomes fractional Langevin equation and its solution 
has been considered by several authors 0, [24Tj26l . |44| . 
The variance and MSD of the position process at long 



3 



time behaves like normal diffusion, that is a 2 ~ t and 
R 2 ~ t as t — > oo. On the other hand, we have <r 2 ~ 



t 



-> 00. 

2 as t 



oo. Note that if (|T0| is replaced 



by the fractional velocity v(i) = D°x(t), < (3 < 1, 
then for x = one has R 2 ~ i 2/3_1 [23. This leads to 
SFD subdiffusion when (3 = 3/4. For t — > 0, one obtains 
i^ 2 t 2 ( Q +^)- 1 . which gives normal diffusion if (3 = 1 — a, 
and ballistic motion ifa< = /3 = 3/4 [11] . 



Case 2. < a < 1, 7(4) = 72 (t) = A 2 r( t 1 J 7) , and 
C(i) = c c t- c , c c > 0, < C < 1 

This is the case of fractional generalized Langevin equa- 
tion without external force and with power law type of 
memory kernel p], d [H . The MSD corresponds to 
this case satisfy the following asymptotic properties. As 
t — > 0, one gets 



if 



t 2 , if a > C/2 

i 2 + 2 "-?, if a < C/2 



(19) 



In the case when a > C/2, the particle undergoes ballistic 
motion for very short times. However, when a > C/2, 
normal diffusion with R 2 ~ i is possible only if C = 1 + a 
or C > 1, which contradicts our assumption that C < 1. 
For t — > 00, 



R 2 



(t^-(, if 7 >C/2 
jlogt, if T = C/2 , 

I constant, if 7 < C/2 



(20) 



which gives subdiffusion with R 2 ~ when 7 = (2C + 
l)/4. Here we note that the large time asymptotic be- 
havior of MSD is independent of the fractional order a of 
the generalized fractional Langevin equation. From (|19[) 
one notices that the particle does not diffuse normally 
for short times, though it can undergo subdiffusion for 
7 > C/2 after sufficiently long time. However, if the 
particle is assumed to undergo ballistic motion at very 
short times, then the generalized Langevin equation in 
this case gives the correct short and long time limits of 
the MSD for SFD. 



III. FRACTIONAL GENERALIZED LANGEVIN 
EQUATION WITH EXTERNAL FORCE 

In this section we consider the generalized Langevin equa- 
tion @ with external force under the following general 
setting: 

Case 3. < a < 1, 7(4) = Ai<5(t) + A 2 rn^) and a ^ °- 
A special case with a = 1, 7 = 1/2 and a = has 
been considered recently in [38| . Solution to the position 
process is then given by (fT4"f with G(t) is given by the 
inverse Laplace transform of 



1 



G{s) = — n r- 

v ; S Q +! + Ai 



s + A 2 st ' 



(21) 



In order to study the asymptotic properties of the solu- 
tion, we expand G(s) in the following series form 



G(s) 



s a+l-1 + \ : 



X 

n=0 



^"7 



s q+1 -t + A 2 



(22) 



One can express the solution in terms of the two parame- 
ter Mittag-Leffler function E a: p(z) by using the following 
Laplace transform relation: 



L [tfE a ,p (Xt a )] = 
where E ai p(z) is defined by [Io| 

E a ,p{z) = 



A' 



n=0 



T(an 



Now define 



G (t) — t a E a _ 1+ i^ a+ i (— A 2 t Q 7+1 ) 
G Q (t)* = t a 1 E a — t +i >a (— A 2 t" 7+1 ) 



(23) 



(24) 



(25a) 
(25b) 



t pui 

du\G* a (t — u\) J du 2 G* (m - u 2 ) ■ 



Jo 

U„-2 



and 
GT(t) 

Jo 

■■■ du n -iG* {u n -2 - Un-l) G* (Un-l) 

Jo 

(26) 

for n > 2. We then have 



Git) = G (t) + V (-Aif / duG Q 

n=0 Jo 



(t~u)G* n (u) (27) 



Some properties of G(t) which are necessary for obtaining 
the asymptotic limits of MSD are given in [A] 

We shall also need the following asymptotic expansion 
of Mittag-Leffler function [l(j to obtain the asymptotic 
properties of the variance and MSD. For z — > 00, 



/v 



(-l) n - 1 z~ n 



0(\z\ 



-1-N\ 



|arg(z)|< (I--)tt (28) 



and for z — > 0, 



(i). Short time limit 
For t — > 0, we have 



O(z), 



(29) 



G(t) 



r 



j2a — i+l 



-Ax: 



Hor(« + l r(2a-7 + 2) T(2a+1) 

(30) 



4 



Since 7 < 1, one gets 



G(t) 



t° 



Ai: 



t-toT(a+l r(2a + l)' 
From dH]) and (|A1|) in one gets 



X Xo — Vol 



Co 



"G(i) + o/ 1 - K G(t) 



t - Ai 



+ a 



T(a + 2) 

a — k-\-1 



T{a-n + 2) Al T(2a-K + 2) 
We therefore obtain 



x~x ~ w t + a- 

t^o r(o;-K + 2) 



(31) 



(32) 



(33) 



(34a) 



and 



Vot, 



a = 0, or a > k 



t->o 



(v + a)t, fl^0,Q = K 
a r(a-l+2) ' «^0,a<K 



(34b) 



Using the following expression for variance (see |B] for its 
derivation) 

a 2 = 2k B T duG(u) - J duG(u)D a G(u) J , (35) 



one gets 



2k B TX 1 



(2a + l)r 2 (a + l) 



(36) 



Here we remark that in the derivation of (|35[) , we have 
made use of the fluctuation-dissipation theorem (see [5J . 
It has been pointed out that the FD theorem fails in the 
presence of external random noise (2lj . However, the ex- 
ternal force in the fractional generalized Langevin equa- 
tion (|9|) is a non-random force, so FD theorem is applica- 
ble in this case up to first order of the external force term 
41 . Note that we also assume that the FD theorem ap- 
plies to the fractional generalized Langevin equation, as 
done in references [28|, [38| . Usual fluctuation-dissipation 
theorem is valid for SFD system provided the fluid or 
gas is elastic and diluted [40[ • In the case of strongly in- 
clastic and dense systems, fluctuation-dissipation formula 
fails and a more general form of fluctuation-dissipation 
relation has to be used. 

By noting that the short time behavior of variance is 
of the order t 2a+1 , therefore for a > 1/2, it is larger than 
(x — x ) given by the square of P4b[) . thus it is only this 
term that dominates the short time limit of MSD. This 
gives 



R 2 ~ { (v + a) 2 t 2 , 

1 V 2 (a-n+2) ' 



a = or a > k 
a^0,a<K 



(37a) 



The short time limits associated with a > 1/2 given by 
(|37a[) implies that the two cases with a > n (zero exter- 
nal force) and a — k (non-zero external force) are both 
ballistic in nature. The case corresponds to non-zero ex- 
ternal force with a < K is non-ballistic, and it leads to 
normal diffusion if n = a + 1/2. 
For a = 1/2, one gets 



R 2 




4fc B TAi \ .2 

7T / ' 

2 , 4k b TXi 



TV 



r%3/2- K )' 



a = or 1/2 > k 
a^0,K = l/2 . 
a/0,l/2<K 

(37b) 



The short time limits given by (|37bj) are quite similar 
to those in (|37a[) . Again, the first two cases with a = 
1/2 > k (zero external force) and a = 1/2 = k (non- 
zero external force) both lead to ballistic motion, while 
the third case with a = 1/2 < k (zero external force) is 
non-ballistic, it gives normal diffusion only if n = 1. 
Finally, for a < 1/2 one obtains 



R 2 



2fc s TA! 



,2 f 



(2a+l)r a (a+l) 

2fc-f 2 



r 2 (a-K+2) ■ 



a = or a > k 
a/0,a < k 



(37c) 



Both cases in (|37c[) can not be ballistic. The first case 
with a > k or zero external force gives superdiffusion, 
hence does not describe SFD. For non-zero external force 
with a < k, the short time limit leads to normal diffusion 
if « = a+1/2, which is the same condition as for the third 
case in (|37a[) and (|37b|) . In other words, we always obtain 
normal diffusion in the short time limit for k = a + 1/2. 

(ii). Long time limit 

For t ^> 1, we have from (l33l) 



X Xo ^ Vo 
t— > 00 



When n > a 



A 2 r( 7 -a + l) 
becomes 



X Xo rsJ Vo 



t n/-a 



From the variance relation (|35[) one gets 
a 2 ~ 2fc s r( f- 

t^oo B \A 2 r( 7 + i) 

^27 — a — 1 



A 2 r( 7 - k+ 1) 

(38) 



(39) 



Al(2 7 -a-l)r( 7 )r( 7 -a) 



2k B T 



V 



t^oo a A 2 r( 7 + i)' 



(40) 



since 7 < a + 1. Using (TT5)) for the MSD and together 
with (|39|) and (|40|) . we obtain the following three cases 
for the long time limit for MSD: 



5 



If 7 < 2a, 



If 7 = 2a, 



R 2 



t^oo A 2 r(7 + i) 



2k R T 



A 2 r( 7 + i) A2r 2 (h + !) 



(41a) 



t 7 . (41b) 



If 7 > 2a, then we have 



If 



t 2( 7 -a) 



X 2 T 2 (7-a + l) 



(41c) 



We remark that the initial velocity can be determined by 
the equipartition principle of kinetic energy v 2 = k B T. 
One may therefore assume that coefficient a for the ex- 
ternal force is proportional to \Jk B T and takes the form 
a = ai\/k B T, a± is a positive constant. Note that the 
fractional order of the Langevin equation or a does not 
appear in the long time limits of MSD for the position 
process in the case of 7 < 2a and the case 7 = 2a, just 
like for Case 1 and Case 2 (without external force) dis- 
cussed in Section [XT] However, when 7 > 2a, the long 
time limit of MSD R 2 given by (|41c[) is dependent on a. 
On the other hand, the exponent of the external force k 
is absent in (T4*Tj) . 

Here we would like take note of the importance of ini- 
tial conditions on SFD which has been pointed out in a 
recent work Q. Here we recall that according to equa- 
tions (|36p and (|37|) . the short time limits for the variance 
and MSD have different time dependence. If now we sup- 
pose the initial condition is x(0) = instead of x(0) = x , 
and we assume the process is centred, that is x = 0, then 
we have R 2 = a 2 . Now the short time limit of MSD has 
the same time dependence as the variance. For the long 
time limits, both the variance and MSD have the same 
time dependence except for 7 > 2a, and they are inde- 
pendent of x(0). A similar observation had been obtained 
in reference [291 ] 

The above discussion shows that it is possible to obtain 
the correct short and long time behaviors of the SFD us- 
ing the generalized fractional Langevin equation ([9]) and 
(|10p under certain specific conditions. In the following 
we want to consider some limiting and special cases. 

Case 3a. Overdamped case: 

Consider the overdamped case such that the Newton ac- 
celeration or ballistic term (for a = 1) or the fractional 
acceleration term (for a / 1) is neglected. Now the sec- 
ond term in the expression for variance (|35|) will be ab- 
sent, and the mean is simply 



(xi(t)} = x , 



such that 



R 2 = a 2 = 2k B T [ duG(u). 
Jo 



(42) 



(43) 



Let us calculate the overdamped case in more detail. 
Now the Laplace transform of G(t) is 



G(s) 



1 



s (l-7)"l 



Ais + A 2 st "V-ff + A' 
where £ = 1/Ai and A = A2/A1, thus 

G(t) = CE7i_ 7 ,i (-At 1 " 7 ) . 

One now gets 

R 2 = a 2 = 2k B T( / dr J Bi_ 7 . 1 (-At 1 " 7 ) 
= 2/c s TC J Bi_ 7 ,2 (-At 1 " 7 ) . 
The asymptotic behaviors of MSD are given by 



(44) 



(45) 



R z 



2k B T(t, 
2k B T( 1 



Ar(i- 7 ) 



t 

t — > 00 



(46) 



(47) 



Therefore for 7 = 1/2, one gets the correct subdiffusion 
behavior corresponds to the SFD. 

Case 3b. < a < 1, j(t) = A 2 p*i- 7 ) ana ^ a 7^ 
This is the case of (|9]) with external force and sin- 
gle term of memory kernel (with Ai = 0). We have 
G(t) = G (t) = t a E a ^ ha+1 (-y a ~ y+ ) from 11253 
as the Green function. Using (|A2[) and (| A3|) in [5] to- 
gether with the asymptotic formulas of Mittag-Lcfflcr 
functions (|2"8|) and (|2"9")) . one gets from equations in ([33)) 
the asymptotic properties of the variance. For t — > 00 



2k B T 



P 



_A 2 r( 7 + l) A|(27-l-a)r( 7 )r(7-a) 

A 2 r(7 + 1) 

where we use the fact that 7 < 1 + a. On the other hand, 
for t -> 



a 2 ~ 2k B T 



(2a - 7 + 2)r(a + l)r(a -7 + 2) 



(49) 



Now consider the mean 



x - x = v Q tE a - 7+1 , 2 (-X 2 t a 7+1 ) 
+ at 



K+1 E a ^ 1+1 ^ K+2 (-A 2 i«- 7+1 ) (50a) 
£ 7 — a tn~ K 



a- 



t^oo A 2 r( 7 -a + l) A 2 r( 7 -K+1)' 

(50b) 

which decays to zero for a > 7 and k > 7 as t — > 00. 
The MSD behaves asymptotically in the same way as 
the variance, i.e. 



R 2 ~ 2k B T 



A 2 r( 7 + i)- 



t — > 00. 



(51) 



6 



One gets from (|50a[) the small time limit 

^a— k+1 

x - x ~ v t + a—, — , t-tQ. (52) 

1 (a — k + 2) 

Since < a < 1 and < k < 1, if a < k, we have 

lOt — K+l 



t -> o. 



r(a-K + 2 

We can now conclude that the MSD takes the form 



(53) 



R 2 ~a 2 



r 2 (a- k + 2) 



2k B T 



A 2 



{2a - 7 + 2)r(a + l)r(a -7 + 2) 



(54) 



The case a = 1/2, 7 = 1/2, k = 1 gives the MSD as 



\2yfc B r 



r 2 (3/2) 



A 2 r(3/2) 







00 



(55) 



which satisfies the behavior of SFD. 

As a manifestation of the non-passing many-particle 
syatem, a realistic physical model of SFD should include 
the long-range inter-particle correlation, and hence the 
particle density term. However, in our model the corre- 
lation function does not have a closed analytic form, we 
do not pursue it since SFD can be characterized by its 
MSD, which can be calculated more easily. The particle 
density and others many-body effect have been absorbed 
into the memory kernel term and the random force, and 
it can be linked to the variance and MSD via the diffu- 
sion coefficient D and SFD mobility F (see for example 
(38|). We shall identify in next section the constants in 
the memory kernel with the diffusion coefficient and SFD 
mobility. The long-range correlation property of SFD in 
our model is manifested in the realization of the single-file 
subdiffusion process as fractional Brownian motion with 
Hurst index H = 1/4. The long-range dependence for 
such a realization can be verified by using the argument 
given in (27j . 

We remark that the mathematical modeling of SFD 
based on its MSD will not be unique. Recall that a 
Gaussian random process is determined (up to a mul- 
tiplicative constant) by its mean and covariance func- 
tion. Two different Gaussian processes can have the same 
MSD or variance but with different covariance functions. 
Thus, any modeling based on fractional Langevin equa- 
tions that give the correct short and long time limits to 
the MSD of SFD can at best be regarded as one of the 
possible models. In order to have a more realistic and 
concrete model of SFD, it is necessary to take into ac- 
count the physical interactions and interplays between 
the boundary conditions and the particle motion. 



IV. NEW CLOSED EXPRESSIONS FOR MSD 
OF SFD 

From the above discussion we see that it is not possible 
to obtain a close analytic expression for the MSD for 
the position process from the solution to © and (|10l) . 
even though we are able to show that both the long time 
and short time asymptotic properties of the MSD agree 
with that for SFD. However, for the overdamped case one 
has the closed expression (|46|) for the MSD. By letting 
7 = 1/2, (Uf)]) gives the correct short and long time limits 
for the MSD of SFD. 

Let us consider this special case in more detail. If we 
let 7 = 1/2, and with appropriate values for f and A, 
we can then show that R 2 = 2k B T(tE 1 / 22 (—At 1 / 2 ) can 
provide the correct description of SFD. By comparing 
this expression with ([3]) one gets in the limit f — > 0, 



1?- 



2k B T(t = l 2 (l - 6) 



f 



(56) 



which requires C = ^ 2 (1 - 0)/(2k B Tr). On the other 
hand, by comparing it with (U) one has for t — > 00, 



R 2 ~ 2k B T(t- y 



r(3/2) 



4fcgTC , 1/2 
A^ 



f 1 — 



(57) 



which requires A = 9^2/t. Therefore, the expression 
R 2 = 2k B T(tE 1/2<2 (-At 1 / 2 ) with C = l 2 {l-6)/{2k B TT) 

and A = 9^2/t provides a new alternative expression for 
the MSD of SFD. 

By noting that both £1/2,2 {t) and Ei/ 2i p(t), (3 > 
have the same short and long time limits (up to a mul- 
tiplicative constant), it is then possible to generalized 
the above closed expression for the MSD to the following 
more general expression 



R 2 — 2k B TQ'tE 1 / 2 f 3 



X't 1 ' 2 



(58) 



for /3 > 0. Again, by comparing (|58|) with <|3j> and Q one 
gets C' = l 2 {l-0)T(p)/{2k B TT) and A' = 6T(/3-±)/V27. 
Expression ([55)1 gives the family of closed form expres- 
sions for SFD, and they can be regarded as alternatives 
to the equation ([58} obtained by Brandani Q. Simula- 
tions of R 2 = 2C t tEi/ 2 2 (—At 1 / 2 ) and the more general 
case R 2 = 2QtE\i 2 ^ (—At 1 / 2 ) with some specific values 
of C, A and /3 are given in Figures Q] and [2] 

The SFD model having transition from normal diffu- 
sion to subdiffusion regime can only be regarded as an 
approximation to the real SFD system. At very short 
times such a system first undergoes ballistic motion with 
MSD ~ t 2 . Ballistic motion occurs before the particle has 
had a chance to collide with anything during the initial 
very short times. There exists the possibility of the direct 
transition from the ballistic regime to the single-file be- 
havior [3 113, EH . Such a tendency becomes more promi- 
nent with increasing particle concentration, and this has 
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been demonstrated by molecular dynamics simulations 
[3 [lH . In order to describe such a situation, one con- 
sider the generalized Langevin equation ((9]) with f(t) = 
and 7(t) = Aa^ 1 / 2 , which gives the MSD as 



R 2 = 2k B Tt 2 E 3/2>3 (-A 2 i 3 / 2 ) 



The detailed derivation of (|59|) is given in [Cj It can be 
shown easily that (|59j) gives the correct long and short 
time limits for a SF system that goes directly from ballis- 
tic motion to the SF regime. Simulation of (|59[) is given 
in Figure [21 We also show that the presence of the three 
regimes, namely the ballistic motion, normal diffusion 
and SFD subdiffusion in Figure [3] based on (1ST 



V. CONCLUDING REMARKS 

We have studied a class of generalized Langevin equa- 
tion with and without external force. It is found that 
with a specific choice of external force, namely one which 
varies as t~ K , < k < 1, then the solution of the general- 
ized Langevin equation © gives a Gaussian process with 
MSD that has the same asymptotic properties as that of 
SFD. 

We remark that an equivalent way to characterized 
SFD is to use the probability distribution function 
W(x, t) which can be obtained as solution to the Fokker- 
Planck equation. Derivation of Fokkcr-Planck equaion 
corresponds to the fractional generalized Langevin equa- 
tion (|9]) is a highly non-trivial problem. Furthermore, it 
may not be fruitful to consider such an equation since we 
only use ((9J to obtain the asymptotic properties of SFD. 
However, it may be interesting to note that W(x, t) at 
the short and long time limits are respectively given by 
solutions to the following Fokker-Planck equations: 



aw(x,t) = D fw(x,t) ^ 



dt 



dx 2 



and 



dW(x,t) F d 2 W(x,t) 



dt 



2y/i dx 2 



t —> oo. 



(60) 



(61) 



(|60p is just the usual diffusion equation and its solu- 
tion gives the probability distribution of normal diffu- 
sion, while (pTTj) is similar to the effective Fokkcr-Planck 
equation for fractional Brownian motion with Hurst in- 
dex H = 1/4 [42|. One may heuristically combine these 
two equations into one: 



dW(x,t) 
dt 



D Vi 



Pit) 



D, 



d 2 W(x,t) 
dx 2 



where 



P(t) 



if t 

1 if t 





oo 



(62a) 



(62b) 



(1621) has the disadvantage that the intermediate values of 
p(t) is not known. One concrete version of Fokkcr-Planck 
equation that gives the short and long time limiting equa- 
tions (fBT)]) and (jrTTj) is the following: 



(59) dW{x,t) 



dt 



/2,1 



2D 



1/2 \ d 2 w(x,t) 

dx 2 



(63) 



which has the required properties since 



E 



and 



E\/2 



1/2,1 



2Do 
F 



2D 



1/2 



1 



t-s-o r(i/2) 



a/2 



-r ( 64 ) 

\/7T 



F 



r(i/2)2&tva %/^D Q t i n 



(65) 



Note that l|6"3")) is one possible effective Fokker-Planck 
equation which gives the correct asymptotic probability 
distributions for SFD. Our remarks on the limitation of 
using MSD to characterize SFD apply to probability dis- 
tribution as well. W{x, t) does not describe SFD uniquely 
even if we know its values for all intermediate times, or 
the correct Fokker-Planck equation for the process. On 
the other hand, if one has the correct fractional Langevin 
equation for SFD at all times, then it describes the pro- 
cess uniquely. 

We also obtain a class of new closed expressions for 
MSD of SFD by considering the solution of the over- 
damped case, thus providing an alternative expression 
to that of Brandani From another special case of 
generalized fractional Langevin equation (0), we derive 
a closed expression describing MSD of SFD with three 
regimes, namely ballistic motion, normal diffusion and 
sub-diffusion. These closed analytic expressions of MSD 
are given in terms of Mittag-LefHer functions. Finally, 
we remark that despite of numerous studies carried out 
on SFD, there are very few results which deal with pos- 
sible concrete realizations of SFD as a specific stochastic 
process. In a recent work, we proposed one such possi- 
ble realization of SFD as the step fractional Brownian 
motion, which has the flexibility of giving the correct de- 
scription to SFD with two and three regimes [29j . 



Appendix A 

Here we give some properties of G(t) which arc required 
for obtaining the asymptotic values of MSD. We have the 
following properties of fractional integrals and derivatives 
ofG(i): 



I l -»G{t) = I 1 ~ fl G (t) 

oo „ t 

+ Y\ (-Ai)™ / duI^G, 

n=l J ° 



{t-u)G* n {u) 
(Ala) 
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and 



D+G(t) =D<t>G (t) 

oo „t 

+ y"(-Ai) n / duD*G {t-u)GT{i 

n=l ^ 

From (|25ap the following can be calculated 



(Alb) 



D a G (t) = £ Q _ 7+M (-A 2 < Q -^ +1 ) (A2a) 

I l -<*G {t) = t£„- 7 +i, a (-A 2 t"- 7+1 ) (A2b) 

J K G {t) = t a K+ i? Q - 7 +i. Q - K +2 (— A2t a 7+ ) • 

(A2c) 

In addition, one needs to consider the following: 



duG (u) = / duu a E a ^ 1+ i, a+ i (-A 2 u Q 7+1 ) 
./o 

= £" + 1 -Ea-7+l,a+2 (~ A 2 £" 7+1 ) ■ (A3) 



Appendix B 

We give the derivation of the expression for variance ([33]) . 
From (fT6]) . the Laplace transform of the convolution 



is given by 



A(t) = / duC(t - u)G(u) 
Jo 



A{s) = C(s)G(s). 



(Bl) 



(B2) 



From the fluctuation dissipation relation we have C (t) = 
k B T^y(t), and express according to definition (|15j) . we 
have 



A = k B T- 



7(s) 



-7(5)5 



1 - 



s a + 7(s) 



(B3) 



Replace the convolution term in (|16j) by the inverse 
Laplace transform of (|B3[) we get ([33)) . 



process is given by 
x(t) = x Q + 



r(i - qJo 



(i-u)~ Q G(w)<iu 



+ / G(t - u)£(u)du, 
Jo 

where G(i) is the inverse Laplace transform of 

1 



G(s) 



s a+l _|_ ^( s ) s 

For 7(t) = A 2 r( f 1 J' 7) , one gets 

G(t) = t Q i? Q _ 7 +i, Q -|-i (— A 2 t Q 7+1 ) . 
The mean and variance of x(t) are respectively 
x = x +v I 1 ~ a G{t) 

and 



(CI) 

(C2) 

(C3) 
(C4) 



({x{t) -5) 2 ) = 2k B T J duG(t)~J duG(t)D a G(u 

(C5) 

Now using the relation R 2 — (x — x ) 2 + <j 2 , one sees that 
the MSD R 2 can have a close form provided the second 
term of the variance cancels out the term (x — x Q ) 2 ■ This 
is possible if a = 1. We thus obtain 

R 2 = 2k B T [ duG{u) = 2k B T [ duuE 2 - 1 a (-A 2 m 2 ~ 7 ) 
Jo Jo 

= 2k B Tt 2 E 2 _^ 3 (-A 2 i 2 " 7 ) . (C6) 
When 7 = 1/2, one gets 



R 2 = 2k B Tt 2 E z/2 ^ (-X 2 t^ 2 



(C7) 



Appendix C 

Consider the fractional generalized Langevin equation @ 
and (p~0|) with f(t) = 0. The solution for the position 
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Long Time 




(a) (b) 
FIG. 1: Short- and long-time simulations for MSD 2tE 1/2 , 2 (-* 1/2 )- 



Short Time 





(a) (b) 

FIG. 2: Short- and long-time simulations for MSD 2t 2 E 3/2 . 3 (-t 3 / 2 ). 




FIG. 3: Plots showing the three regimes of SFD according to MSD R 2 = 2t 2 £ , 3/2i3 (-< 3 / 2 ). 



